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Abstract. Lattice simulations of QCD iiave produced precise estimates 
for the masses of the lowest-lying hadrons which show excellent agree- 
ment with experiment. By contrast, lattice results for the vector and 
axial vector form factors of the nucleon show significant deviations from 
their experimental determination. We present results from our ongo- 
ing project to compute a variety of form factors with control over all 
systematic uncertainties. In the case of the pion electromagnetic form 
factor we employ partially twisted boundary conditions to extract the 
pion charge radius directly from the linear slope of the form factor near 
vanishing momentum transfer. In the nucleon sector we focus specifi- 
cally on the possible contamination from contributions of higher excited 
states. We argue that summed correlation functions offer the possibility 
of eliminating this source of systematic error. As an illustration of the 
method we discuss our results for the axial charge, gx, of the nucleon. 



1 Introduction 

Simulations of QCD on a space-time lattice have recently succeeded in producing 
reliable results for several phenomenologically important quantities. This enormous 
progress resulted mainly from improvements of numerical techniques, which led to a 
significant acceleration of simulation algorithms. A variety of hadronic observables can 
now be computed with fully controlled statistical and systematic errors. A well-known 
example is the spectrum of the lowest-lying mesons and baryons. Following years of 
efforts by many different collaborations [l|2)3)4|5|6|7|8|9)10j . it has now been firmly 
established that QCD accounts for the experimentally observed spectrum within the 
quoted uncertainties of the lattice calculation, which are at the level of a few percent. 
This underlines once more that QCD is the correct theory of the strong interaction, 
also in the low-energy regime. 

Surely, the role of lattice QCD is not restricted to the verification of known re- 
sults. For instance, the masses of the light quarks are not directly accessible by ex- 
periment but can be predicted by lattice calculations. Recent years have witnessed 
the publication of a wealth of lattice results for the strange quark mass, as well as 
the isospin-averaged light quark mass. The availability of accurate predictions for a 
number of phenomenologically relevant observables has prompted the foundation of 
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Fig. 1. Compilation of results for the mass of the strange quark from the FLAG report [llj . 
The top panel shows results from simulations with dynamical up, down and strange quarks 
{N{ = 2 + 1), while the middle panel refers to QCD with two light flavours. In the bottom 
panel (blue points) results from non-lattice determinations are plotted. Red symbols denote 
results that have been published in refereed journals. Results from preprints or conference 
proceedings are represented by green points. The global estimate for QCD with Ai'f = 2 + 1 
flavours is represented by the grey band, while the corresponding estimate for Nf = 2 flavours 
is bounded by the dashed lines. 



the FLAG Working Group which tries to form global averages for quark masses, me- 
son form factors and decay constants, and other quantities, very much in the spirit of 
the Particle Data Group. Figure[T] shows a compilation of results for the mass of the 
strange quark from the FLAG review . A remarkable feature is the impressive level 
of consistency, despite the strong variation of systematic effects among the different 
determinations. FLAG'S analysis of all published lattice data leads to the "global" 
estimates of JJ[i 

mfS(2GeV) = 94±3MeV, mjf|(2 GeV) = 3.43 ± 0.11 MeV, (1) 

where rriud — + nij). Here the masses are quoted in the MS-scheme of dimen- 

sional regularisation, at the commonly used reference scale of2GeV. This level of 
precision exceeds that of the values quoted in the current edition of the Particle Data 
Book [12] by an order of magnitude. 

Among the other quantities discussed in the FLAG review are form factors for 
Ki3 decays and the ratio /k/Ztt of kaon and pion decay constants. When combined 
with the experimentally measured branching fractions for the respective leptonic and 
semi-leptonic decays, the lattice estimates can be used to test the unitarity of the 
first row of the CKM matrix. As discussed in detail in [TT] , the test can be further 
strengthened by including the constraint of the direct determination of \Vud\ from 
nuclear /3-decay. In this way, first-row unitarity is confirmed at the permille level, 
using only lattice results and experimental data as input. 
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Fig. 2. Recent compilation of lattice results for the nucleon axial charge, plotted versus the 
pion mass squared (from ref. [13). The various sets of symbols refer to data from different 
collaborations. The left-most point denotes the value from the Particle Data Book. 



In spite of these successes, one finds examples for which the agreement between lat- 
tice calculations and experiment is less satisfactory. Observables that describe struc- 
tural properties of the nucleon fall into this category, as was pointed out in several 
recent reviews on the subject [13ll4ll5j . For instance, the experimentally observed 
dependence of the (isovector) electric and magnetic form factors of the nucleon on 
the squared momentum transfer is not reproduced. Moreover, lattice calculations 
of the nucleon axial charge, ^a, lie typically 10 — 15 % below the experimental value, 
as signified by the summary plot in Fig.[2j What is even more disturbing is the ab- 
sence of a clear trend in the lattice data which would indicate that the gap becomes 
narrower as the pion mass is decreased towards its physical value. Since one would 
hesitate to conclude that QCD has been falsified on the basis of these observations, 
the contradiction can only be resolved if one accepts that systematic effects in lattice 
calculations of these observables are not fully controlled. 

The goal of our project is the calculation of a variety of nucleon and meson form 
factors with complete control over all systematic uncertainties. As we shall see, the 
application of novel techniques which allow for a more reliable identification of the 
ground state is a crucial ingredient for our task. 

The outline for the remainder of this article is as follows: In section [2] we review 
the basic "lattice technology", including a discussion of the main systematic effects 
inherent in any lattice calculation. Section [3] describes our on-going project aimed 
at producing a benchmark calculation of the pion electromagnetic form factor. In 
section |4] we discuss our determination of the nucleon form factors and the axial 
charge. Finally, section [5] contains our conclusions. 

2 Lattice technology 

Every lattice simulation proceeds by performing a stochastic calculation of observables 
via Monte Carlo integration. If the gluon field is represented by the link variables 
U^{x), the expectation value of some observable fi is given by 
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where dC/^(x) denotes an integration over the group manifold of the gauge group 
SU(3). Each flavour / = m, d, s,... contributes a factor of the quark determinant, 
where is a suitable representation of the massless Dirac operator. It is important 
to bear in mind that the discretisation is not unique: Standard discretisations such 
as Wilson or staggered fermions have been known for a long time. At the end of 
the 1990s, alternative lattice transcriptions of the quark part of the action based on 
the Ginsparg- Wilson relation were shown to preserve chiral symmetry at non- 
zero lattice spacing [17118] , and particular realisations such as domain wall |19l20j 
or overlap fermions |21) have been applied in practical simulations. More recently, 
several different implementations of so-called "minimally doubled" fermions have been 
studied {22 23 24 25 26 27 , in the hope of finding discretisations which have good 
chiral properties but avoid the large inherent numerical costs of domain wall or overlap 
quarks. 

A common feature of all discretisations is the strong growth in the computer 
time required for generating statistically independent configurations of gauge fields 
as the light quark masses are tuned towards the physical values of the up- and down- 
quark masses [55]. During the past 10 years this algorithmic problem has been greatly 
ameliorated owing to several significant technical improvements. These include hi- 
erarchical integration schemes [29130131] , mass preconditioning '35^ , domain decom- 
position methods [3D] , deflation techniques [33] and suitably optimised combinations 
thereof [34]. The value of the pion mass, i.e. the lightest mass in the pseudoscalar 
meson channel, serves as a measure for how deeply a particular simulation has pene- 
trated into the chiral regime. As a result of the recent algorithmic improvements one 
can now routinely access pion masses as low as 200 MeV, while in 2001 that figure 
stood at a heavy 500 MeV. In a few cases, contact with the physical pion mass has 
already been made [35l36j . 

We will now discuss the main systematic effects in lattice calculations. For any 
non-zero value of the lattice spacing a, observables receive corrections of order a^, 
where the integer p depends on the chosen discretisation. Thus, the value of p deter- 
mines the rate of convergence towards the continuum limit. In practical simulations 
the continuum limit is taken by computing observables for several values of a before 
performing an extrapolation to a = 0. The relatively high numerical cost of simula- 
tions in the chiral regime implies that results at the physical pion mass are obtained 
via chiral extrapolations. Although the chiral behaviour of many observables can be 
constrained by Chiral Perturbation Theory (ChPT), such extrapolations are still a 
major source of uncertainty. This is particularly true for quantities which describe 
structural properties of the nucleon. 

Observables computed on the lattice are also affected by the finite spatial vol- 
ume. Empirically, one finds that uncertainties arising from the finite box size L are 
subdominant, provided that L is at least as large as 2.5 — 3fm and that the pion 
mass in units of L satisfies rriTrL ^3 — 4. Whether this is true for any given observ- 
able must be established on a case-by-case basis. In fact, there are suspicions that the 
above bounds are not sufficient to guarantee small finite-size effects for many baryonic 
quantities [14j . It should also be noted that finite- volume effects can be computed an- 
alytically in ChPT, which offers the possibility of applying finite-volume corrections 
to the final lattice estimates. 

An important systematic effect arises from the possibility that masses and matrix 
elements of ground state hadrons are contaminated by contributions from excited 
states. Physical observables are usually extracted from correlation functions of local 
composite fields. If Ohad denotes an interpolating operator for a particular hadron. 
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Fig. 3. Two spatial profiles of smeared quark fields obtained using the Jacobi algorithm |37] 
in combination with HYP-smeared link variables f3B j . The width of the profile is controlled 
by the number of iterations and the value of the mass parameter in the kernel of the smearing 
function. 



then its two-point correlation function is given by 



x,y n 

(3) 

where the sum over n arises due to the insertion of a complete set of eigenstates of 
the Hamiltonian. The quantity Wn{p) denotes the spectral weight of the nth state. 
The operator Ohad projects on all states that are characterised by the same quantum 
numbers. The last relation shows that the ground state dominates the correlation 
function for large Euclidean time separations, (?/o ~ xq) 3> 0, which allows for the 
determination of the ground state energy Ei (p) from the exponential fall-off. If, how- 
ever, the statistical fluctuations grow rapidly with increasing (?/o — a;o) it may happen 
that the signal is lost before the contributions from excited states are sufficiently sup- 
pressed. As a result, one incurs an uncontrolled distortion of ground state properties. 
A widely used procedure which is designed to enhance the projection onto the ground 
state (i.e. which increases its weight wi{p) in the spectral representation) is called 
"smearing". Common smearing algorithms apply a kernel function F{x,y; U) to the 
quark field at point y, in order to approximate the spatial profile of the hadron's wave 
function. Two examples of smeared sources are shown in Fig.|3] 

Our simulations have been performed on the high-performance cluster "Wilson" 
operated by the Institute of Nuclear Physics at Mainz University. We use N{ = 2 
flavours of 0{a) improved Wilson fermions as our discretisation of the quark action. 
The non-perturbative estimates for the improvement coefficient Csw which multiplies 
the Sheikholeslami-Wohlert term [B^ were taken from ref . [10] . Monte Carlo ensembles 
were generated for three different values of the lattice spacing, i.e. a w 0.08, 0.07 and 
0.05 fm and for a range of pion masses from TOtt ~ 250 MeV to 700 MeV. Lattice 
volumes were chosen sufficiently large so as to satisfy m^L > 4 on all ensembles. A 
compilation of simulation parameters is shown in Table[T] The Monte Carlo ensembles 
are being generated as part of the CLS ("Coordinated Lattice Simulations") projectF] 



^ https : //twiki . cern . ch/ twiki/bin/view/CLS/WebHome 



6 



Will be inserted by the editor 



Table 1. Simulation parameters and approximate values for the lattice scale and pion 
masses. The preliminary results presented in this review are based on the ensembles labelled 
"A", "E", "F" and "N". 
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Fig. 4. Assignment of twist angles for the three-point correlation function of the vector 
current. 



3 The pion electromagnetic form factor 

The electromagnetic form factor, defined by 

(7r+(p^)||Tl7^w- p7^rf|7r+(p,)) = {pf + Vi)t. (4) 

where q — pf—pi is the momentum transfer, encodes the distribution of electric charge 
inside the pion. Of particular interest is the charge radius, (rj), which is derived from 
the pion form factor at vanishing momentum transfer, i.e. 



6^ "'^ ' ' dq2 



(5) 



Lattice calculations of mesonic matrix elements are technically simpler than the cor- 
responding quantities for the nucleon. Furthermore, the pion electromagnetic form 
factor receives no contributions from quark-disconnected diagrams, whose evaluation 
typically suffers from large statistical fluctuations. This opens the possibility to per- 
form a precision test of lattice QCD, by comparing lattice estimates for (r^) to the 
experimentally determined value. However, owing to the finite spatial volume the 
accessible range of momentum transfers is severely constrained, which presents a 
major obstacle for precise lattice determinations of (rj). 

The solution to this problem, which by now has become a standard method, is 
to employ partially twisted boundary conditions |41|42j . To this end one imposes 
periodicity on the quark fields up to a general phase, i.e. 

^(x + LeW) =e^^''V(^), fc = l,2,3, (6) 
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Fig. 5. Left: pion form factor computed for a range of pion masses compared to the results 
of [45146) . Right: data points from the inset in the top left-hand corner. 



where e^'^^ denotes a unit vector in the k^^ spatial direction, and 6'^*''-' is the corre- 
sponding phase angle. A non-zero value of the latter modifies the accessible values of 
the spatial momentum according to 



p = n 



27r 

T 



e 

L' 



(7) 



Quark propagators computed for different values of the twist angles can be combined 
to form the three-point correlation functions from which the pion form factor can be 
determined (see Fig.|4|. The twist angles for the initial and final state pions are then 
given by [JS] 



Oi — 6i — 6 



3- 



Of —0-2 



9. 



so that the expression for the squared momentum transfer becomes 



Pf 



L 



L 



(8) 



(9) 



Thus, by an appropriate choice of twist angles one can tune to any desired value. In 
our simulations we have chosen 0^ , Of so as to achieve a particularly fine momentum 
resolution near — Q. 

So far we have assumed that boundary conditions are identical for sea and va- 
lence quarks. It is, however, customary to apply partially twisted boundary condi- 
tions, where the twist is applied only to the quark fields in the valence sector. This 
has the advantage that the generation of Monte Carlo ensembles must be performed 
only once (e.g. for zero twist), while the choice of twist angle and, in turn, the mo- 
mentum transfer can be optimised for a particular observable. The modification of 
the boundary conditions can lead to finite-size effects associated with the breaking 
of flavour symmetries. However, it was shown in ref. [44j that such finite-size effects 
are exponentially suppressed in processes without final-state interactions. Hence, for 
the electromagnetic interaction between a photon and a pion finite-size effects are 
expected to be small. 

In Fig. [5] we show our results for the pion form factor computed on the ensembles 
N3, N4, N5 and F6. There are two main observations: First, there is a clear trend 
towards a steeper fall-off with q^ = —Q^ as the pion mass decreases from about 
600 MeV on N3 to about 290 MeV on F6. Secondly, by our choice of twist angles 
we were able to produce a very dense set of points near = 0. This allows us to 
extract the pion's charge radius in an accurate and model-independent fashion, by 
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Fig. 6. The squared pion charge radius (in units of ro) extracted from the linear slope of 
the form factor in an interval [0, (gro)^], plotted versus the interval length. The meaning of 
the labels is given in Table[l] 
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Fig. 7. The squared pion charge radius as a function of the squared pion mass. The black 
open and solid symbols are taken from ref. |45| . The value from the Particle Data Book is 
indicated by the pink star. All dimensionful quantities are expressed in units of the hadronic 
radius ro. To locate the positions of the pion masses in physical units we have set ro = 0.5 fm. 



determining the linear slope of over a narrow interval, starting at 5^ = 0. 

For the following discussion we express all dimensionful quantities in units of the 
hadronic radius tq |47l48j . In Fig.[6]the values of (rj) /rg are plotted versus the length 
of the interval in (qro)"^ over which the slope was determined. Obviously one would 
like to choose this interval as small as possible. The figure shows that the statistical 
accuracy in the determination of (r^)/rQ is still very good in the immediate vicinity 
of vanishing momentum transfer. The fact that the resulting estimates of (r^) are 
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practically constant implies that terms of 0{q^) in the chiral expansion of /^r are 
quite small. 

Our preliminary results for the squared charge radii were determined from the 
linear slope over the interval —0.15 < (qro)'^ < 0. Using tq = 0.5 fm this corresponds 
to \q\ ^ 150 MeV. FigureM shows the chiral behaviour of the charge radius. By com- 
paring the solid red and olue points we conclude that our data are accurate enough 
to exhibit a sensitivity to lattice artefacts. Overall, though, the chiral trend compares 
favourably with the phenomenological value of (rj), shown as the pink symbol in the 
plot. A more systematic investigation of lattice artefacts, as well as a more detailed 
study of the g^-dependence of the form factor is left for future work. 



4 Nucleon form factors and axial charge 

The pion form factor discussed in the previous section can be considered a warm- 
up exercise for the technically more difficult case of extracting the corresponding 
quantities in the nucleon sector. The well-known Dirac and Pauli form factors, Fi and 
F2 , are related to the matrix element of the electromagnetic vector current between 
nucleon initial and final states via 



{N{p',s')\V^,ix)\N{p,s))^u{p',s') j^Fiiq') + a,,,^F2iq') u{p,s) 



(10) 



where \N{j), s)) denotes the initial state of a nucleon with momentum p and spin s. 
The electric and magnetic form factors Ge and Gm are derived from 



G^{q^) ^ F,{q^) 



(2mN) 



r^2(<z'), GM{q^)^F,{q^)+F2{q^). 



(11) 



Similarly, the matrix element of the axial current A°(a;) = is pa- 

rameterised in terms of the form factors Ga and Gp, i.e. 



{N{p\ s') \Al{x)\N{p, s)) = u{p', s'y^T'^ [l,l5GA{q') + |^Gp(q2)] uip, s). (12) 

The axial charge, qa, is defined as the axial form factor at vanishing momentum 
transfer, gA = hm Ga('?^). Nucleon matrix elements are determined by computing 

the corresponding three-point correlation functions. After performing the Wick con- 
tractions one can express the correlators in terms of quark propagators. In general, 
this gives rise to the quark-connected and quark-disconnected diagrams depicted in 
Fig.|8] Since the disconnected contribution is statistically very noisy, it is mostly ne- 
glected in lattice calculations, which is the approach we have adopted as well during 
the first stage of this project. 

Unlike the case of the pion form factor there have been only a few attempts to 
apply partially twisted boundary conditions to calculations of nucleon form factors. 
This more cautious approach is motivated by the observation that the induced finite- 
volume effects can become sizeable for small twist angles 43] ■ 

It was already mentioned in the introduction that lattice calculations have so far 
failed to reproduce the experimentally observed dependence of electromagnetic form 
factors on the momentum transfer. This, in turn, implies that the associated charge 
radii are not consistent cither. Moreover, lattice calculations tend to underestimate 
the axial charge qa by 10— 15 %. Uncontrolled systematic effects are held responsible 
for this. With only a few exceptions, nucleon form factors have mainly been computed 
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Fig. 8. Quark-connected and -disconnected contributions to three-point correlation func- 
tions of the vector and axial-vector currents between nucleon states. 



over a very limited range of lattice spacings |50l51l52l53l54l55l56l57l58l59l60j . It is 



thus conceivable that the agreement between lattice calculations and experiment is 
improved when performing a systematic continuum extrapolation to a = 0. Also, 
nucleon form factors could be more sensitive to finite- volume effects, compared to their 
mesonic counterparts "14' . Another potentially very important source of systematic 
error is the chiral extrapolation. The quark masses used in current simulations are still 
relatively large and leave a long extrapolation to the physical value of the (isospin- 
averaged) light quark mass. In particular, the region around the physical pion mass is 
not well constrained by the data, and the functional form for the extrapolation may 
be too crude to yield reliable results. 

What makes the calculation of nucleon form factors a much more difficult task 
compared to the corresponding mesonic quantities, is the relatively high inherent level 
of statistical noise in baryonic correlation functions. While the signal-to-noise ratio 
of pseudoscalar meson correlators can be shown to remain constant as the Euclidean 
time separation is increased, one finds that the noise increases exponentially in the 
baryonic sector [61162163] . The high level of statistical fluctuations is the source of 
another potentially very dangerous systematic effect, namely the contamination of 
ground state properties by contributions from excited states, which have not died out 
in the region of Euclidean times where statistical errors are still quite small. Unless 
one has interpolating operators at one's disposal which maximise the spectral weight 
of the ground state (see eq. ([s])) one either risks the distortion of results from excited 
state contributions or has to accept large statistical errors, provided that the signal 
is not lost entirely. The commonly used smearing methods to enhance the spectral 
weight of the ground state may be insufficient to guarantee reliable results. They must 
be combined with more sophisticated techniques such as the generalised eigenvalue 
problem [ 64165166] . or, as in our project, summed operator insertions [67168] . 

In order to illustrate our approach we restrict ourselves to the discussion of the 
nucleon axial charge. For several reasons, gA is an ideal benchmark quantity for lattice 
calculations of structural properties of the nucleon: (1) is derived from a matrix 
element of a simple quark bilinear which contains no derivatives, (2) initial and final 
nucleon states are at rest, and (3) its definition as an iso- vector quantity implies that 
contributions from quark-disconnected diagrams are absent. 

Nucleon form factors are usually extracted from suitably chosen ratios of three- 
and two-point functions, such as 



Rr{q;t,ts) 



C[{q,t,ts) f C2{q,ts-t)C2{0,t)C2{0,ts) 
C2(0, ts) i C2(0, ts - t) C2iq, t) C2(q, ts) 



> 1/2 

I , (13) 
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where F ~ V, A characterises the current which is inserted at time t. The nucleon 
is created from the vacuum at time zero and annihilated at time tg > t. For gj^, the 
momentum transfer q — p' — p is zero. Assuming that the axial current is correctly 
renormalised and that all kinematical factors are properly taken into account the 
ratio i?A gives direct access to ^a, i-e. 

Rit, t,) = RA{q = 0, t, tg) = 5A + 0(c-^*) + 0(e-^(*=-*)). (14) 

Here, A denotes the mass gap between the ground state nucleon and its first excita- 
tion. In QCD with dynamical quarks A is expected to be equal to 2171^^, which implies 



that the corrections to in eq. ( 14 1 can be rather sizeable, unless large values of t 
and ts are considered. However, since the statistical errors grow exponentially with 
time separation, it is difficult to optimise the choice of t and t^, in order to keep both 
the systematic errors due to excited state contamination and the statistical errors un- 
der control. Most published results for ^a were obtained by fitting the ratio R{t,ts) 
to a constant in t, for « 1 fm. Below we present evidence that this procedure may 
be insufficient to ensure that the resulting estimates of ^a are free from excited state 
contributions. 

Let C^{t,tg) denote the three-point correlation function of the axial current for 
q — 0. For a fixed value of tg we define the ratio p{t, tg) as 

'^'^'^^ SMry ^-«'2a,...,t.. (15) 

Obviously, the deviation of p(t;ts) from unity at a particular value of t indicates 
the presence of excited state contributions. If the three-point functions are computed 
using smeared sources (and perhaps also sinks) the deviation of p{t,ts) from one can 
be regarded as a measure of the effectiveness of the smearing procedure. Since the 
axial charge is usually determined by fitting R{t,ts) to a constant for t around ts/2, 
the parameters of the smearing procedure must be tuned such that excited state 
contributions are eliminated inside the fit range. 

In Fig. [9] the ratio p{t, tg) computed using smeared-local (SL) correlator^ for tg = 
1.1 fm is plotted against t. While p is mostly compatible with one for small values of 
t, there are large and significant deviations for t > is/2. One concludes that source 
smearing is unable to remove excited state contamination in the interval < t < tg. 
The problem is further highlighted by comparing p{t,ts) to the effective mass of the 
nucleon: The bottom panel in Fig.|9] shows clearly that the asymptotic behaviour of 
the nucleon two-point function only sets in at timeslices t for which the deviation of 
p from unity becomes significant. We can draw the conclusion that source smearing 
alone cannot guarantee the absence of contamination from excited states in ratios such 
as R(t, ts). Our findings suggest that there is a mismatch in the asymptotic behaviour 
of three- and two-point functions which enter the ratio R{t,ts). As a consequence, 
the appearance of a plateau in R{t, tg) for t, tg 1 fm must considered to be mostly 
accidental. One might expect that the situation improves if correlators are smeared 
both at the source, t = 0, and sink, t = tg. This is currently under investigation. 

The use of summed operator insertions [ST] offers a handle for eliminating excited 
state contributions. The key observation is that the corrections from excited states in 



eq. (14) can be parametrically reduced. Restricting the discussion once more to the 
case or the axial charge, one defines the summed ratio S{tg) via 

S{tg):=Y,R{t,tg). (16) 



^ The notation "SL" is used to describe a correlator which is smeared at the source, t = 0, 
only. 
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Fig. 9. The ratio p{t,ts) (top panel) and the effective mass of the nucleon (bottom panel), 
computed on a lattice of size 32'' • 64 (L = 2.2 fm) and a pion mass of 415 MeV. Two- and 
three-point functions were computed using Jacobi-smeared sources. The parameters in the 
smearing procedure were chosen to produce a smearing radius of ~ 0.5 fm. 



Its asymptotic behaviour is given by 

S{t,) = c + gAts + Oit,e-^'^), (17) 

where the (divergent) constant c includes contributions from contact terms. Since 
tg > t the corrections to S{ts) are in general more strongly suppressed than for 
R{t, tg). The reduction of excited state contributions in S{ts) comes at a price, though, 
since the summed ratio must be computed for several different values of ts, in order to 
extract from the linear slope. As a further comment we add that the method can 
be straightforwardly extended to cases where the initial and final states have different 
momenta. The general expression is given as 

Sriq;Q -.^Y. ^riq;t,ts) = K + Miq^% + OiUe-'^'^) + 0{t,e-^'''), (18) 
t=o 

where M(q^) is the matrix element of interest, and A, A' denote the energy gaps in 
the initial and final states, respectively. 

In our simulations we have computed the summed correlator ratios for ~ 0.7 — 
1.1 fm, using both the vector and axial vector currents. As an illustration how the 



method works, we compare in Fig. 10 the t-dependence of the ratio R{t, tg) computed 
for three different values of ts against the result extracted from the linear slope of 
the summed correlator. The latter is shown as the purple band in the figure. It is 
seen that — 15a « 1.1 fm is just sufficient to produce a plateau which agrees with 
the result from the summed insertion within statistical errors. Nonetheless, it is clear 
that the summed correlator produces a larger value for ^a- 

The case for using summed correlators is even more compelling as the chiral limit 
is approached. Since A 2m.n. decreases for lighter quark masses, the correction 



terms in eq. (14) become larger and may spoil the expected chiral behaviour of g^. 
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Fig. 10. The ratio R(t,ts) of the iso-vector axial charge computed for three different values 
of ts a,t P = 5.3 (a ~ 0.07 fm). The purple band denotes the result extracted from the slope 
of the summed correlator. 



1.4 



1.2 



O.f 



Prelirhinary 



a = O.OB fm 

• a = 0.07 fm 

• a = 0.05 fm 



0.5 



1 1.5 



Fig. 11. The axial charge qa determined using summed correlators plotted versus the pion 
mass squared. The magenta star represents the PDG value of gA ~ 1.2695(29). 



The size of the corrections in the standard approach also depends on the spatial extent 
L of the lattice, since the overlap of a local interpolating field with a multi-particle 
state is suppressed by powers of the volume. In the conventional approach it is then 
difficult to disentangle finite-volume effects from excited states contaminations. A 
slight drawback of the method can be read off from Fig.[TO] Excited state contributions 
to the summed correlator are reduced at the expense of incurring larger statistical 
errors. 
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The current status of our gA determination is summarised in Fig. |ll| All data 
points were determined via the slope of the summed correlators which results in 
larger statistical errors compared to conventional calculations using similar statistics. 
Nonetheless, it is clear that these preliminary results are in good agreement with 
the experimental value for pion masses ^ 300 MeV. These findings differ from 
what is usually observed by other collaborations. Clearly, further studies including 
data produced at smaller pion masses and additional sets of correlation functions 
which are smeared both at the source and the sink are required before any definite 
conclusions can be drawn. Also, we will perform a systematic investigation of the 
infiuence of lattice artefacts, by analysing the results obtained at three values of the 
lattice spacing. 



5 Summary and conclusions 

Despite the recent successes of lattice QCD it is clear that quantities describing struc- 
tural properties of hadrons are still afflicted with one or several sources of systematic 
errors. In this work we have outlined and applied a number of technical improvements. 
In particular, summed insertions have proven a valuable tool to suppress excited state 
contamination in ratios of correlation functions. Moreover, partially twisted boundary 
conditions lead to a much enhanced momentum resolution in calculations of the pion 
form factor, which greatly facilitates the extraction of the pion's charge radius with 
reduced model dependence. Finally, our large lattice volumes and fine lattice spacings 
ensure that the corresponding systematic effects are under good control. 

Owing to its computational simplicity, the pion form factor is an ideal testbed for 
the more complicated case of nucleon form factors. Clearly, more work, including the 
calculation of correlators which are smeared both at the source and sink, is required 
before any definite statement about the computed g^-dependence and its comparison 
with experimental data can be made. The axial charge of the nucleon, ^a, is of partic- 
ular importance, since its determination does not involve the technically challenging 
calculation of disconnected diagrams, and because it is defined at a single value of 
. Both of these features make it an ideal reference quantity for future benchmark 
calculations of structural properties of the nucleon. 
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